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SUMMARY 


A procedure for the numerical solution of the complete isothermal elasto- 
hydrodynamic lubrication problem for rectangular contacts is outlined. This 
procedure calls for the simultaneous solution of the elasticity and Reynolds 
equations. In the elasticity analysis the conjunction is divided into equal 
rectangular areas. It is assumed that a uniform pressure is applied over each 
area. In the numerical analysis of the Reynolds equation the parameter $ = 
Qhp/2 -j s introduced in order to help the relaxation process. The analysis 
couples the elasticity and Reynolds equations, going from the inlet to the 
outlet without making any assumptions other than neglecting side leakage. 

By using the procedures outlined in the analysis the influence of the 
dimensionless speed U, load W, and materials G parameters on minimum film 
thickness is investigated. Ten cases are used to generate the minimum- 
film-thickness relationship. 


H . = 3.07 U°’ 71 G°* 57 W 0,11 

min 

The most dominant exponent occurs in association with the speed parameter; the 
exponent on the load parameter is very small and negative. The materials pa- 
rameter also carries a significant exponent, although the range of the param- 
eter in engineering situations is limited. The five dimensionless speed 
parameter values used in obtaining the preceding equation are varied over a 
range six times the lowest speed value. The four dimensionless load values 
are varied over a range 1.8 times the lowest load value. Conditions corre- 
sponding to the use of solid materials of bronze and steel and lubricants of 
paraffinic and naphthenic mineral oils are considered in obtaining the expo- 
nent in the dimensionless materials parameter. 

Plots are presented that indicate in detail the pressure distribution, 
film shape, and flow within the contact. The characteristic pressure spike is 
clearly in evidence as is the parallel film shape through the central portion 
of the contact. Minimum film thickness occurs near the outlet of the contact. 


INTRODUCTION 

The recognition and understanding of elastohydrodynamic lubrication rep- 
resents one of the major developments in the field of tribology in the twen- 
tieth century. The revelation of a previously unsuspected regime of lubrica- 
tion not only explained the remarkable physical action responsible for the 
effective lubrication of many nonconformal machine elements like gears and 
rolling-element bearings, but also brought order to the complete spectrum of 
lubrication regimes, ranging from boundary to hydrodynamic. 

Historically one of the first to study the effect of elastic distortion 
in highly loaded contacts was Meldahl (ref. 1). He examined the effect of 
high load on film shape and pressure distribution for a constant-viscosity 
lubricant. The first attempt to analyze both elastic and viscous effects in 
elastohydrodynamical ly lubricated contacts was done by Grubin and Vinogradova 
(ref. 2), who managed to incorporate both the effects of elastic deformation 
and the viscosity-pressure characteristics of the lubricant in the inlet anal- 
ysis of hydrodynamic lubrication of nonconformal contacts. Their work dealt 
mostly with a line contact, and it was assumed that the shape of the elastic- 
ally deformed solids in highly loaded lubricated contacts was the same as the 
shape produced in dry (Hertzian) contacts. This assumption facilitated the 



solution of the Reynolds equation in the inlet region of the contact and en- 
abled the separation of the solids in the central region of the contact to be 
determined with commendable accuracy. 

Dowson and Higginson (ref. 3) obtained an empirical formula for the iso- 
thermal-line-contact elastohydrodynamic lubrication problem. This formula 
showed the effect of speed, load, and material properties on minimum film 
thickness and was based on their theoretical solutions. In the procedure they 
adopted, the computed film shape was compared with the shape of the elastic- 
ally deformed solids, and then the pressure curve was modified to improve the 
agreement between the two shapes. These calculations were performed on hand- 
operated desk calculating machines. Jacobson (ref. 4) solved the elastohydro- 
dynamic lubrication problem for a spherical contact under pure rolling condi- 
tions and a non-Newtonian lubricant with a limiting shear strength. The shear 
strength of the lubricant limited the pressure gradients and shear stresses in 
the oil. Therefore no pressure spikes were seen in the theoretical solutions. 

Hamrock and Dowson (ref. 5) were able to successfully obtain a theoreti- 
cal approach to coupling the elasticity equation with the Reynolds equation 
for elliptical contacts such as those normally found in gears and ball bear- 
ings. Hamrock and Dowson 1 s work on elastohydrodynamic lubrication (ref. 6) 
considered the complete spectrum of contact geometries (ranging from point to 
line contacts), materials (hard and soft), and lubricant availability (fully 
flooded or starved conditions). 

To obtain a better understanding of the failure mechanism in machine ele- 
ments, the next generation of elastohydrodynamic lubrication analysis should 
incorporate such effects as 

(1) Surface roughness effects 

(2) Non-Newtonian effects 

(3) Temperature effects 

The foundation studies that produced the complete elastohydrodynamic lubrica- 
tion solutions for rectangular contacts given in the present report will be 
used in further studies incorporating these effects. To incorporate these 
effects in an initial study, rectangular contact analysis, often referred to 
as "line" or "one dimensional" contact analysis, should be used instead of 
elliptical contact analysis because of the added complexity of considering 
these effects. 

The analysis in the present report couples the elasticity and Reynolds 
equations, going from the inlet to the outlet without making any assumptions 
other than neglecting side leakage. This analysis may thus be used as the 
foundation for the more complicated analysis incorporating surface roughness 
effects, non-Newtonian effects, and temperature effects. 

In the results the influence of dimensionless speed, load, and materials 
parameters on minimum film thickness was investigated for a contact fully 
immersed in lubricant (i.e., fully flooded). The dimensionless speed and load 
parameters were varied over a range of 11 and 2-1/2 times, respectively. Con- 
ditions equivalent to using solid materials of bronze and steel and lubricants 
of paraffinic and naphthenic mineral oils were considered in obtaining the 
exponent on the dimensionless materials parameter. Ten cases were used in 
obtaining the fully flooded film thickness formula. A fully flooded condition 
is said to exist when the inlet distance of the contact ceases to influence 
the minimum film thickness in any significant way. The inlet distance of the 
contact is defined as the distance from the center of the contact to the edge 
of the computing area. Besides the film thickness calculations that were 
made, calculation of the force components, shear forces, coefficient of fric- 
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tion, and center of pressure were also performed. A simple formula is pro- 
vided that described the location of the center of pressure relative to the 
inlet of the contact and as a function of the dimensionless load and speed 
parameters for steel surfaces. 

Computer plots are presented that indicate in detail the pressure spike 
and minimum film thickness. Plots are also shown of the parallel shape of the 
reduced pressure, the Poiseuille term, and the mass flow per unit length with- 
in the contact. 


SYMBOLS 

A 

B 

b 

b 

C 

Ci , . . . , Cg 

D 

E 

E 1 

F 

f 

G 

H 

^min 

H min 


constant defined in eq. (27) 

bib 

semiwidth of Hertzian contact, 2R V2 W/tt, m 
b/n, m 

constant defined in eq. (28) 
constants 

defined by eq. (19) 

modulus of elasticity, N/m2 

effective elastic modulus, 

2/{[(l - vl)/E a ] + [(1 - vg)/E b ]}, N/m2 

dimensionless shear force 

shear force per unit length, N/m 

dimensionless materials parameter, aE ' 

dimensionless film thickness, h/R 

dimensionless minimum film thickness, h m -j n /R 

dimensionless minimum film thickness obtained from 
least-squares fit of data 


(H . ) 


min' 


D 


(H . ) 
min 7 


DH 


dimensionless minimum film thickness obtained from Dowson 
(ref. 7) 

dimensionless minimum film thickness obtained from Dowson 
and Higginson (ref. 3). 

dimensionless constant defined in eq. (23) 

film thickness, m 
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minimum film thickness, m 
dimensionless mass flow per unit length 
mass flow per unit length, kg/(s m) 
constant defined in eg. (29) 
constant defined in eq. (30) 
number of nodes in semi-axis of contact 
dimensionless pressure, p/E 1 
pressure, N/m^ 

dimensionless reduced pressure, q/E' 

reduced pressure, N/m^ 

effective radius in x-direction, m 

curvature radius, m 

geometrical separation, m 

defined in eq. (67) 

dimensionless speed parameter, nou/E'R 

surface velocity in direction of motion, (u a + u^)/2, m/s 

(^min - Hnviri) 100/H m -j n 

[(H . ) - H. 1100 /TT . 

L min q min- mm 

dimensionless load parameter, w z /E'R 
load, N/m 

dimensionless coordinate, x/b 

location of dimensionless center of pressure, x C p/b 

coordinate in direction of motion, m 

location of center of pressure, m 

coordinate in direction of film thickness, m 

pressure-vi scosity coefficient of lubricant, m^/N 

angle defined in eq. (43) 

elastic deformation, m 
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Subscripts: 

a 

b 

x 

z 


coefficient of determination 

absolute viscosity at gage pressure, (N s)/m2 

dimensionless viscosity, n/no 

viscosity at atmospheric pressure, (N s)/m^ 

coefficient of friction 

Poisson's ratio 

lubricant density, kg/m^ 

dimensionless density, p/pq 

density at atmospheric pressure, kg/m^ 

QH 3 / 2 

solid a 
solid b 

coordinate in direction of motion 
coordinate in direction of film thickness 


THEORY 

Reynolds Equation 

The general approach to the numerical solution of the one-dimensional 
rectangular or line-contact problem covered in this report is similar to the 
method used by Hamrock and Oowson (ref. 5) in solving the two-dimensional 
elliptical contact problem in elastohydrodynamic lubrication. The Reynolds 
equation for one-dimensional flow where side leakage is neglected can be 
written as 



where u = (u a + %) /2 is the mean surface velocity or the entraining 
velocity in the x-direction. 

Letting 


X = x/b, p = p/pQ, n = n/n q, H = h/R, and P = p/E' (2) 
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where 



equation (1) can be rewritten in dimensionless form as 


(3) 

(4) 



(5) 


( 6 ) 


is the dimensionless speed parameter, 

W 7 

W = TR (7) 

is the dimensionless load parameter, and w z is the load per unit width. 

Figure 1 shows the radius of the rollers used in defining equation (3). 

It is assumed that convex surfaces, as shown in figure 1, exhibit positive 
curvature and concave surfaces, negative curvature. Therefore, if the center 
of the curvature lies within the solid, the radius of curvature is positive; 
if the center of curvature lies outside the solid, the radius of curvature is 
neg ati ve. 


r b 


r a 


(a) Two undeformed rollers. 

(b) Equivalent cylindrical solid near plane. 

Figure 1. - Rollers and equivalent roller. 
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The viscosity of a fluid may be associated with the resistance to flow, 
with the resistance arising from intermolecular forces and internal friction 
as the molecules move past each other. Because of the large pressure varia- 
tion in the lubricant films the viscosity of the lubricant is not constant in 
elastohydrodynamical ly lubricated contacts. Barus (ref. 8) proposed the 
following formula for the isothermal viscosity-pressure dependence of liquids: 

n = noe a P (8) 


In dimensionless form this equation can be written as 



(9) 


where 
G = aE 1 

is the dimensionless materials parameter. 

Substituting equation (9) into equation (5) gives 


( 10 ) 


4r(»H 3 e~ GP £) - 2W < U > 

Grubin and Vinogradova (ref. 2) were the first to introduce the writing of the 
pressure and viscosity in terms of a reduced pressure as 


0 = | (i - e" GP ) 

and 

dQ -GP dP 

dx “ e dx 


(12) 


(13) 


Note that, as P » °°, Q * 1/G. Substituting equation (13) into equation (11) 
gives 



(14) 


Density 

For a comparable change in pressure the density change is small as com- 
pared with the viscosity change. However, very high pressures exist in 
el astohydrodynamic films, and the liquid can no longer be considered as an 
incompressible medium. From Dowson and Higginson (ref. 9) the dimensionless 
density for mineral oil can be written as 


p 


-e- « l + 

p o 


0.6 E'P 
1 + 1.7 E'P 


( 15 ) 


where E' 


is given in gigapascals. 



Film Shape 

The film shape can be written simply as 
h(x) = h 0 + S(x) + «(x) (16) 

where 


h 0 constant 

S(x) separation due to geometry of undeformed solids 
6(x) elastic deformation 


The separation due to the geometry of the two undeformed rollers shown in 
figure 1(a) can be described by an equivalent cylindrical solid near a plane, 
as shown in figure 1(b). The geometrical requirement is that the separation 
of the two rollers in the initial and equivalent situations should be the same 
at equal values of x. Therefore using the parabolic approximation we can 
write the separation due to the undeformed geometry of the two rollers as 


S(x) 



(17) 


Figure 2 shows a rectangular area of uniform pressure width 2b._ From 
Timoshenko and Goodier (ref. 10) the elastic deformation at a point x on the 
surface of a semi-infinite solid subjected to a pressure p at the point 
xi can be written as 


«(x) 




p In (x - Xj) 


2 



(18) 


Since the pressure is assumed to be uniform over the rectangular area, the 
pressure can be put in front of the integral. The integration of equation 
(18) then results in the following: 


6 ( x) = | PD (19) 

where 



Figure 2. - Surface deformation of semi-infinite body subjected 
to uniform pressure over rectangular area. 
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D = b[(X - B) ln(X - B) 2 - (X + B) ln(X + B) 2 + 4B(1 - In b)] (20) 

and 

b semiwidth of Hertzian contact 

B b/b 

¥ b/n 

n number of nodes within semiwidth of Hertzian contact 

Now the_term 6(7) in equation (19) represents the elastic deformation at a 
point x due to a rectangular area of uniform pressure p and width 2b. If 
the contact is divided into a number of equal rectangular areas, the total 
deformation at a point x due to the contributions of the various rectangular 
areas of uniform pressure in the contact can be evaluated numerically. The 
total elastic deformation caused by the rectangular areas of uniform pressure 
within a contact can be written as 




( 21 ) 


where 


j = | k - i | + 1 (22) 

Therefore, substituting equations (17) and (21) into equation (16) while 
writing the film thickness in dimensionless form gives 



(23) 


Phi ($) Solution 

Having defined the density and film thickness, we can return to the solu- 
tion of the Reynolds equation. The dimensionless reduced pressure Q given 
in equation (14) plotted as a function of X exhibits a very localized pres- 
sure field with high values of dQ/dX and d^Q/dX^. Such a condition 
with high gradients is not welcomed when performing numerical analysis by 
relaxation methods. Therefore, as discovered by Hamrock and Dowson (ref. 5), 
to produce a much gentler curve, a parameter $ is introduced, where 

♦ = qh 3/2 , e- GP l ( 24) 

The dimensionless reduced pressure Q is smaller at large values of film 
thickness H and large when the film thickness is small. The 4> substitu- 
tion also has the advantage that it eliminates all terms containing deriva- 
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tives of products of H and Q or H and $. Therefore, from equation 

(24) while expanding terms within equation (14), we get 

» 3/2 -sr (» &) - r- 1(» H 1/2 $)- 24U < 25 > 

The first step is then to replace equation (25) by a finite difference 
approximation. The relaxation method relies on the fact that a function can 
be represented with sufficient accuracy over a small range by a quadratic 
expression. With standard finite central-difference representation, equation 

(25) can be written as 


. - L.$. - M. = 0 
l i+l l i-l i i i 


(26) 


where 


C i ~ p i+l + 3p i-l 
L i = 4(p i+l + p i-l } 


(27) 

(28) 


, K5[p i + i V^Ti^m - 4H i + H i-i> + p i-i - 4H i + 3H i-i)^ 


H 3 ' 1 2 

i 


(29) 


M. = 
1 


24U VF (l) (o i+l H i+l - 


H 3 ' 2 


(30) 


Figure 3 shows the uniform placement of the nodes within the contact. This 
nodal structure was used in all of the results. The number of nodes within a 
semicontact width was 120 throughout all calculations. 

The following boundary conditions have been adopted: 

(1) At the inlet and outlet the pressure was set to zero. This implies 

that Q and $ are also zero at these positions. 


-Hertzian zone- 

Center of 
contact 


Outlet 


-180 Nodes- 


Inlet 


120 Nodes - 


480 Nodes - 
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— 660 Uniform nodes 

Figure 3. - Nodal structure for numerical calculations. 




Figure 4. - Flow chart of computational procedures used in 
elastohydrodynamic lubrication studies. 


(2) At the cavitation boundary 

p --ay=° (3i) 

Eauation (26) represents a system of simultaneous equations that was solved by 
the standard Gauss-Seidel iterative method. A flow chart shown in figure 4 
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describes the computational procedures used in the elastohydrodynamic lubrica- 
tion studies. 


Mass Flow Rate Per Unit Length 

The mass flow rate per unit length for elastohydrodynamically lubricated 
contacts can be written as 

k = Up h - (32) 


Making use of equation (2) while rearranging terms gives 



Making use of equation (13) allows this equation to be rewritten as 



The second term on the right side of equation (34) is the Poiseuille or pres- 
sure term. Writing the reduced pressure gradient in equation (34) in a cen- 
tral difference form and rearranging the terms, we get 



This expression enables the reduced pressure to be written in terms of flow, 
film shape, and density at the preceding location. 


Force Components 

Figure 5 shows the force components acting on the two solids along with 
the film geometry in a portion of a concentrated contact. Conventionally only 
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the z-components of the normal forces acting on the solids (w az and W£, z ) 
are considered. However, it was felt that the tangential x-components (w ax 
and W{j X ), shear forces (f a and f^), coefficient of friction u, and 
center of pressure x C p should also be expressed and quantitative values 
obtained for each of these expressions. The normal z-component of the force 
per unit length acting on the solids can be written as 

w z = w az = wt, z = /p dx (36) 

Using the dimensionalization given in equation (2) gives this equation as 

M ■ rsr ■ ? / p dx ■ ! (/ p dX ) 2 < 37 > 

In equation (37) W is usually referred to as the dimensionless load param- 
eter. 

The tangential force component w ax for solid a is zero. The tangen- 
tial force component per unit length for solid b is not zero and can be 
written as 

V = - / p dh = - /p-a£ dx (38) 

Using integration by parts gives 

w bx = " [ph:, i + / h 7bT dx (39) 

where i and o refer to inlet and outlet edge of the computing zone, re- 
spectively. 

However, the pressure at the edge of the inlet and at the edge of the 
outlet is zero. Using the dimensionalization given in equation (2) gives 

“bx ■ ftt * / H ^ dx < 40 > 

x 


The resulting force components per unit length can be written as 



Shear Forces 

The shear force per unit length acting on solid a (shown in fig. 5) can 
be written as 




From Hamrock and Dowson (ref. 6, p. 141, eq. 5.30) we can write 

du (2z - h) dp n ^ u a u b^ 
n dz “ 2 dx h 


( 44 ) 


(45) 


Substituting this equation into (44) gives 



Making use of equation (39) while writing equation (46) in dimensionless terms 
gives 


F 


a 





) 



(47) 


The shear force per unit length acting on solid b can be written as 


f b */("&),.„ dx 


Making use of equation (45) gives equation (48) as 


bx 


4(*i - uj 


j- u /\ 3 b § § 

>" r ■<«.*•»> 


Vf/i 


dx 


Note that for equilibrium to be satisfied the following must be true: 

F - F + W, =0 
a b bx 

W - W, = 0 

a z bz 

The coefficient of friction is written as 


u = - 


-F. + W. 
b bx 


(48) 


(49) 


(50) 

(51) 


( 52 ) 


Center of Pressure 

A calculation very useful in traction studies is the location of the 
center of pressure. The appropriate equation is 
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( 53 ) 


X CP -Tr/ px dx 

Writing this in dimensionless form gives 




PX dX 


(54) 


The location of the center of pressure indicates the position on which the 
resulting force is acting. The fact that the resulting force is not acting 
through the center of the roller creates a rolling resistance in the form of a 
moment. This has a significant effect in the evaluation of the resulting 
forces and power loss in traction devices and other machine elements. 


RESULTS 

Dimensionless Grouping 

From the variables of the numerical analysis the following dimensionless 
groups can be defined: 

Dimensionless film thickness: 

H = £ (55) 

Dimensionless load parameter: 


where w^ is the load per unit length. 
Dimensionless speed parameter: 


Dimensionless materials parameter: 

G = ctE ' (58) 

The dimensionless film thickness for a rectangular contact can thus be 
written as a function of the other three parameters: 


H = f(W, U, G) (59) 

The most important practical aspect of elastohydrodynamic lubrication theory 
is the determination of the minimum film thickness within the contact. There- 
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(a) Dimensionless pressure and film thickness. 

(b) Dimensionless reduced pressure and $. 

(c) Dimensionless flow and Poiseuille term. 


Figure 6. - Representative variations of dimensionless 
pressure, film thickness, and reduced pressure, 4>, 
dimensionless flow, and the Poiseuille term on x-axis 
for dimensionless speed parameter U of 1x10"^ 
dimensionless load parameter W of 1. 638x10” , and 
dimensionless materials parameter G of 5000. 
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fore, in the fully flooded results to be presented, the dimensionless param- 
eters (W, U, and G) will be varied and the effect on the minimum film 
thickness will be studied. 

Representative variations of dimensionless pressure, film thickness, 
reduced pressure, *, dimensionless flow, and the Poiseuille term on the 
X-axis are shown in figure 6. Figure 6(a) shows the variation of dimension- 
less pressure and film thickness on the X-axis for U = 1.0xl0-H, W = 
1.638x10-5, and g = 5000. The inlet region is to the left and the outlet 
is to the right in this and each of the remaining figures to be presented. 

The Hertzian pressure is also shown in this figure. The characteristic pres- 
sure spike is clearly evident in this figure as is the parallel film shape 
through the central part of the contact, with a minimum film occurring near 
the outlet of the contact. 

Variation of dimensionless reduced pressure Q and <t> on the X-axis 
for U = 1.0xl0“H, W = 1.638x10“5, and G = 5000 is shown in figure 
6(b). Recall that the reduced pressure is defined in equation (12) and $ in 
equation (24). Figure 6(b) shows that the reduced pressure Q is constant 
within the contact and that $ is constant in part of the contact outlet. 

Figure 6(c) shows the variation of dimensionless flow and the Poiseuille 
term on the X-axis for U = 1.0xl0-H, W = 1.638x10“^, and G = 5000. 

Equation (34) defines the mass flow rate per unit length, and the Poiseuille 
term is the reduced pressure gradient term of equation (34). In figure 6(c) 
the flow is constant throughout the contact. Great care was taken to assure 
that this was true for all the results to be presented. Slight adjustments in 
the pressure profile were necessary in the inlet region to assure that the 
flow was constant in that region. The Poiseuille term approaches 1 at the 
inlet and zero in the Hertzian contact region except for having negative 
values from the pressure spike to the outlet. 


Influence on Load 


Chanqes in the dimensionless load parameter W can be achieved while 
keeping the other parameters constant by changing only the normal applied load 
per unit length w z in equation (56). The values at which the remaining 
parameters U and G were held constant during the calculations were 


U = 1. OxlO-H 
G = 5000 


(60) 


Four values of the dimensionless load parameter W and the corresponding 
values of minimum film thickness H m -j n obtained from the el astohydrodynamic 
lubrication theory developed earlier are shown in table I. These four pairs 
of data were used to determine an empirical relationship between the dimen- 
sionless load and the minimum film thickness. 




(61) 


By applying a least-squares power fit to the four pairs of data [(W-j, 
H m1 - n .,•), i = 1,...,4], the values of Ci and C2 were found to be 
C} =’6.33xl0“6 and C 2 = -0.1056 = -0.11. Therefore the influence of 
load on minimum film thickness can be written as 
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cc 


( 62 ) 


min 


w -°.n 


In addition to the least-squares fit, a coefficient of determination e 
was obtained. The value of e reflects the fit of the data to the resulting 
equation: unity representing a perfect fit, and zero the worst possible fit. 
The coefficient of determination e for these results was 0.9856, which is 
excellent. 

Figure 7 shows the variation of dimensionless pressure and film thickness 
on the X-axis for two values of dimensionless load W (1. 638x10“^ and 
3.0x10"5). The values of the dimensionless speed and materials parameters 
were held fixed as described by equation (60). Figure 7 shows that although 
the pressure distribution changes considerably, the minimum film thickness 
is only slightly affected. This illustrates the slight effect of the dimen- 
sionless load parameter W on minimum film thicknesses, as described in 
equation (62). 


Influence of Speed 


If the surface velocity in the x-direction is changed, the dimensionless 
speed parameter U is modified as shown in equation (57), but the other dimen- 
sionless parameters (W and G) remain constant. The values at which these di- 
mensionless parameters were held constant in the calculations performed to 
determine the influence of speed on film thickness were 


W = 2.048x10-5] 
G = 5000 


(63) 


Values of the dimensionless speed parameter U and the corresponding 
minimum film thickness H mi - n as obtained from the elastohydrodynamic lubri- 
cation of rectangular contacts developed earlier in the report are presented in 


5x10-3 



Figure 7. - Variation of dimensionless pressure and film 
thickness on X-axis for two values of dimensionless 
load. 
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table II. Calculations were performed for five values of dimensionless speed 
parameter covering nearly an order of magnitude. The solutions enabled the 
relationship between minimum film thickness and the speed parameter to be 
written in the form 


H . 


rm n 



(64) 


By applying a least-squares power fit to the five pairs of data [ ( U n - , 

H m -j n i = 1, . .., 5], the values of C 3 and C 4 were found to be 
C 3 =’1179. 2 and C 4 = 0.70640 = 0.71. The coefficient of determination 
e for these results was excellent at 0.9992. Therefore the influence of 
speed on minimum film thickness can be written as 


,0.71 


mi n 


(65) 


Figure 8 shows the variation of dimensionless pressure and film thickness 
on the X-axis for two values of dimensionless speed U (0.5xl0 - H and 
3.0x10“^). The values of the dimensionless load and materials parameters 
were held fixed as described by equation (63). This figure shows that the 
pressure at any location in the inlet region rises in going from the lower 
speed to the higher speed. This result is consistent with the elastohydro- 
dynamic lubrication theory of elliptical contact given by Hamrock and Dowson 
(ref. 5). Note that the pressure spike moves more toward the outlet for the 
lower speed than for the higher speed. A typical el astohydrodynamic film 
shape with an essentially parallel section in the central region is also shown 
in figure 8 . There is a considerable change in film thickness as the dimen- 
sionless speed is changed, as indicated by equation (65). This illustrates 
most clearly the dominant effect of the dimensionless speed parameter U on 
the minimum film thickness in el astohydrodynamically lubricated contacts. 


X 


o 


l.OxlO' 4 4xl0“ 3 



Figure 8. - Variation of dimensionless pressure and film 
thickness on X-axis for two values of dimensionless 
speed. 
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Influence of Material Properties 

A study of the influence of the dimensionless materials parameter G on 
minimum film thickness has to be approached with caution since in practice it 
is not possible to change the physical properties of the materials, and hence 
the value of G, without influencing the other dimensionless parameters con- 
sidered earlier. Equations (56) to (58) show that as either the materials of 
the solids (as expressed in E 1 ) or the lubricant (as expressed in no and 
a) are varied, not only does the materials parameter G change, but so do the 
dimensionless speed U and load W parameters. 

The results obtained from calculations performed for three values of the 
dimensionless materials parameter are summarized in table III. A general form 
of these results, showing how the minimum film thickness is a function of the 
dimensionless materials parameter, is written as 



where 


T = 


H . 
mm 


( 66 ) 


(67) 


In equation (67) the exponents have been rounded off to two significant 
figures so that any error would be absorbed in C 5 , given in equation ( 66 ). 

By applying a least-squares power fit to the three pairs of data, the values 
of C 5 and Cg were found to be C 5 = 3.12 and C 5 = 0.5670 - 0.57. 

The coefficient of determination for these results was 0.9969, which is excel- 
lent. Therefore the effect of the dimensionless materials parameter on minimum 
film thickness can be written with adequate accuracy as 



G°.57 


( 68 ) 


Minimum-Film-Thickness Formula 

The proportionality equations (62), (65), and ( 68 ) have established how 
the minimum film thickness varies with the load, speed, and materials param- 
eters, respectively. This enables a composite dimensionless minimum-f i lmthick- 
ness formula for a fully flooded, isothermal elastohydrodynamic rectangular 
contact to be written as 


H . = 3.07 U°- 71 G°- 57 W 0,11 

mm 


(69) 


In equation (69) the constant 3.07 is different from C 5 = 3.12 mentioned 
earlier to account for the rounding off of the materi al s-parameter exponent. 
In dimensional terms this equation is written as 


h min “ 3 * 07 


0.57 d 0.4, x0.71 

a R ( un Q ) 

/£, \0.03 w 0.11 
z 


(70) 
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From this equation we find that the minimum film thickness depends inversely 
on the effective elastic modulus E' and load per unit length w z . Both 
have small exponents, indicating that the minimum film thickness h mi - n is 
only slightly affected by the effective elastic modulus and load per unit 
length. In contrast to these effects from equation (70) we find that the film 
thickness depends directly on the pressure-viscosity coefficient of the lubri- 
cant a, the geometry R, the surface velocity in the direction of motion u, 
and the viscosity at atmospheric pressure nQ. From the values of the expo- 
nents on these parameters (a, R, u, and no) it is clear that they have a 
dominating effect on the minimum film thickness. 

It is interesting to compare equation (69) with earlier derived minimum- 
film-thickness formulas. Dowson and Higginson (ref. 9) obtained the following 
expression: 


(Tf . ) = 1.6 U°- 7 G°' 6 W °* 13 

mm DH 


( 71 ) 


It was found that this equation produced a positive exponent on the effective 
elastic modulus, which is contrary to physical intuition. This equation was 
therefore revised by Dowson (ref. 7) as 


/TT \ O rr . .0 . 7„0 . 54. .— 0 . 1 3 

Trn'n q ‘ 2 - 65 U G W 


(72) 


The powers of U, G, and W in equations (69) and (72) are quite similar con- 
sidering the different numerical procedures on which they are based. 

The Hamrock and Dowson (ref. 11) ellipical contact minimum-film-thickness 
formula for very long elliptical contacts, where the ellipticity parameter is 
large and the elliptical contact approaches a rectangular contact, is 
expressed as 


(Tf . ) 

min 2 


HD 


3.63 U°' 68 G°’ 49 W 0,073 


(73) 


It is interesting to compare equation (73) with equation (69). The powers of 
U, G, and W are again seen to be quite similar. 

Table IV gives the 10 cases used in obtaining equation (69). In this 
table H m -j n corresponds to the minimum film thickness obtained from the 

el astohydrodynami c lubrication rectangular contact theory developed earlier in 
this report and Tf m -j n is the minimum film thickness obtained from equation 
(69). The percentage difference between these two values is expressed by 
Vi, which is defined as 


'TT 


V 1 " 


min 
TT 


min 1 


100 


min 


(74) 


In table IV the values 
less minimum film thickness 
equation (72) is also shown 


of Vj are within ±2 percent. The dimension- 
obtained by Dowson (ref. 7) and expressed in 
in table IV as (Tf mip ) . The percentage 
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difference between this film thickness and that of H m -j n is expressed by 
V 2 , where 



(TT . ) - F . 

v min' m n 


min 


D 


min 


F . 

min 


100 


( 75 ) 


In table IV the values of V 2 are within +4 to +8 percent, meaning that the 
Dowson (ref. 7) formula produces 4 to 8 percent larger film thickness than 
that obtained from the present analysis. 

Figure 9 compares the pressure profile as obtained from Dowson and 
Higginson (ref. 9) with the present results for U = 1.0x10” , W = 3.0x10”^, 

and G = 5000. The Hertzian pressure is also shown in this figure. The 
Dowson and Higginson (ref. 9) profile is exactly equivalent to the Hertzian 
pressure for most of the contact region, but the present results are somewhat 
lower than the Hertzian pressure. The pressure spike is higher in magnitude 
in the present results and located farther away from the exit of the contact 
than in the results obtained by Dowson and Higginson (ref. 9). 

Table V shows the values of the load components, the coefficient of fric- 
tion, and the location of the center of pressure for the 10 cases presented in 
table IV. The values of the dimensionless load, speed, and materials param- 
eters corresponding to a particular case can be obtained from table IV. Table 
V shows that the tangential force components are three orders of magnitude 
less than the normal force components. Also the coefficient of friction 
decreases with increasing load and increases with increasing speed. Table V 
gives values of the location of the center of pressure for the 10 cases 
evaluated. In all the cases the center of pressure is in front of the center 
of the Hertzian contact. An approximate formula for the location of the 
center of pressure is a function of the dimensionless load and speed and is 
given as 



Figure 9, - Comparison of pressure profile of present results with that 
of Dowson and Higginson (ref. 9). 


(76) 


cp 


= -0.05 U°' 44 W 1,13 


where X C p is measured from the center of the Hertzian contact. This 
formula was obtained by using the data in table V and is valid only for steel 
surfaces. The value obtained from the approximate expression of the location 
of the center of pressure (eq. (76)) is shown in the last column of table V. 
The maximum error in calculating the location of the center from equation (76) 
as compared with the exact solution is 0.89 percent of the Hertzian half 
width, which is excellent. Using equation (76) for bronze surfaces gives an 
error of up to 3.54 percent of the Hertzian half width. This larger error is 
due to having only two cases of bronze material in the formulation of the 
expression. 


CONCLUDING REMARKS 

A procedure for the numerical solution of the complete isothermal elasto- 
hydrodynamic lubrication problem for rectangular contacts has been outlined. 
This procedure calls for the simultaneous solution of the elasticity and 
Reynolds equations. In the elasticity analysis the conjunction was divided 
into equal rectangular areas. It was assumed that a uniform pressure was 
applied over each element. In the numerical analysis of the Reynolds equation 
the parameter $ = QH 2 ' 2 was introduced in order to help the relaxation 
process. The analysis coupled the elasticity and Reynolds equations, going 
from the inlet to the outlet without making any assumptions other than 
neglecting side leakage. 

By using the procedures outlined in the analysis the influence of the 
dimensionless speed U, load W, and materials G parameters on minimum 
film thickness was investigated. Ten cases were used to generate the mini- 
mum-film-thickness relationship 

m- o n7 n0.71p0.57 -0.il 
H min = 3 ' 07 U G W 

The most dominant exponent occurred in association with the speed parameter; 
the exponent on the load parameter was very small and negative. The materials 
parameter also carried a significant exponent, although the range of the 
parameter in engineering situations is limited. The five dimensionless speed 
parameter values used in obtaining the preceding equation were varied over a 
range six times the lowest speed value. The four dimensionless load values 
were varied over a range 1.8 times the lowest load value. Conditions corre- 
sponding to the use of solid materials of bronze and steel and lubricants of 
paraffinic and naphthenic mineral oils were considered in obtaining the 
exponent in the dimensionless materials parameter. 

Plots were presented that indicate in detail the pressure distribution, 
film shape, and flow within the contact. The characteristic pressure spike 
was clearly in evidence as was the parallel film shape through the central 
portion of the contact. Minimum film thickness occurred near the outlet of 
the contact. 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio, August 5, 1982 
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TABLE I. - EFFECT OF DIMENSIONLESS LOAD 


TABLE II. - EFFECT OF DIMENSIONLESS SPEED 


PARAMETER ON MINIMUM FILM THICKNESS 


Force per 
unit length, 

w 7> 

N/m 

Dimensionless 

load, 

W 

Dimensionless 
film thickness, 

^min 

40 000.0 
50 000.0 
60 000.0 
73 249.3 

1.6382xlO' 5 

2.0478 

2.4573 

3.0000 

20.327x 10" 6 
19.711 
19.396 
19.055 


PARAMETE R ON MINIMUM FILM THICKNESS 


Surface 

velocity, 

u, 

m/s 

Dimensionless 

speed, 

U 

Dimensionless 
film thickness, 
Hmin 

0.297040 

0.5xl0 -11 

I2.357xl0‘ 6 

.415856 

.7 

15.482 

.594080 

1.0 

19.711 

1.188160 

2.0 

33.364 

1 .782240 

3.0 

43.029 


TABLE III. - EFFECT OF DIMENSIONLESS MATERIALS PARAMETER ON 
MINIMUM FILM THI CKNESS 


Solid 

material 

Lubricant 

Dimensionless 

materials 

parameter, 

G 

Dimensionless 

speed 

parameter, 

U 

Dimensionless 

load 

parameter, 

W 

Minimum film 
thickness from 
EHL rectangular 
theory, 

H min 

H . 
mm 

u 0 - 7 V 0>11 

Bronze 

Paraffinic 

2553.7 

1 .9579x10'! 1 

4. 0094X10- 5 

20.156xl0- 6 

265.20 

Bronze 

Naphthenic 

3591.1 

5.5975 

4.0094 

52.502 

327.67 

Steel 

Paraffinic 

5000.0 

1 .0000 

2.0478 

19.711 

388.11 
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TABLE IV 


- EFFECT OF DIMENSIONLESS LOAD, SPEED, AND MATERIALS PARAMETERS ON MINIMUM FILM THICKNESS 


Dimensionless 
load parameter, 
W 

Dimensionless 
speed parameter, 
U 

Dimensionless 

materials 

parameter, 

G 

Dimensionless 

film thickness 

obtained from 

EHL theory, 

H . 
min 

Dimensionless 

film thickness 

obtained from 

least-squares 

_fit, 

H . 
min 

Percentage 

difference 

between H . 

min 

and H . , 
mm 

V 1 

Dimensionless 

film thickness 

obtained from 

Dowson (ref. 7), 

(H . ) 
mm D 

Percentage 
difference 
between H 

mm 

and (l W o - 
V 2 

1.6382xl0" 5 

l.OOOOxlO' 11 

5000.0 

20.327x10“® 

20.5 lOx 10 -6 

0.900 

22.020xl0‘ 6 

7.362 

2.0478 

1 



19.711 

20.013 

1.532 

21.391 

6.886 

2.4573 

1 

* 



19.396 

19.616 

1.134 

20.890 

6.495 

3.0000 




19.055 

19.189 

.703 

20.355 

6.076 

2.0478 

.50000 



12.357 

12.235 

-.987 

13.167 

7.617 

2.0478 

.70000 



15.482 

15.536 

.349 

16.664 

7.261 

2.0478 

2.0000 



33.364 

32.737 

-1.879 

34.749 

6.146 

2.0478 

3.0000 

1 

f 

43.029 

43.658 

1.462 

46.154 

5.717 

4.0094 

1.9579 

2553.7 

20.156 

20.420 

1.310 

21.826 

6.885 

4.0094 

5.5975 

3591.1 

52.502 

52.283 

-.417 

54.735 

4.690 



TABLE V. - VALUES OF LOAD COMPONENTS, COEFFICIENT OF FRICTION, AND 
LOCATION OF CENTER OF PRESSURE FOR 10 CASES EVALUATED 


Dimensionless 

D imensionless 

Coefficient 

Dimensionless 

Dimensionless 

normal load. 

tangential load. 

of friction. 

center of 

center of 

W 

W bx 

)J 

pressure, 

X cp 

pressure 
obtained from 
least-squares 
fit, 

J 

cp 

1 .6382xl0' 5 

2.0221x10"® 

6. 1716xl0' 4 

-0.1912 

-0.1881 

2.0478 

2.2939 

5.6008 

-.1552 

-.1463 

2.4573 

2.4071 

4.8979 

-.1240 

-.1191 

3.0000 

2.5448 

4.2414 

-.0971 

-.0951 

2.0478 

1.5388 

3.7572 

-.1041 

-.1080 

2.0478 

1.7647 

4.3087 

-.1194 

-.1251 

2.0478 

2.7985 

6.8330 

-.1893 

-.1982 

2.0478 

3.4243 

8.3609 

-.2317 

-.2367 

4.0094 

4.1764 

5.2082 

-.1031 

-.0921 

4.0094 

7.3268 

9.1369 

-.1810 

-.1456 
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